Nonlinear Damping in Nanomechanical Beam 

Oscillator 
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Abstract — We investigate the impact of nonlinear damping on 
the dynamics of a nanomechanical doubly clamped beam. The 
beam is driven into nonlinear regime and the response is mea- 
sured by a displacement detector. For data analysis we introduce 
a nonlinear damping term to Duffing equation. The experiment 
shows conclusively that accounting for nonlinear damping effects 
is needed for correct modeling of the nanomechanical resonators 
under study. 

Index Terms — Mechanical damping, nonlinear, bistability, 
NEMS, multiple scales. 



I. Introduction 

THE field of micro-machining is forcing a profound redef- 
inition of the nature and attributes of electronic devices. 
This technology allows fabrication of a variety of on-chip 
fully integrated sensors and actuators with a rapidly growing 
range of applications. In many cases it is highly desirable to 
shrink the size of mechanical elements down to the nano- 
scale [1], [2]. This allows enhancing the speed of operation 
by increasing the frequencies of mechanical resonances and 
improving their sensitivity as sensors. Moreover, as devices 
become smaller, their power consumption goes down and their 
cost can be significantly lowered. Some key applications of 
nanoelectromechanical systems (NEMS) technology include 
magnetic resonance force microscopy (MRFM) [3], [4] and 
mass-sensing [5], Further miniaturization is also motivated 
by the quest for mesoscopic quantum effects in mechanical 
systems [6], [7], [8]. 

A key property of systems based on mechanical oscillators 
is the rate of damping. For example, in many cases the 
sensitivity of NEMS sensors is limited by thermal fluctua- 
tion which is related to damping via the fluctuation dissi- 
pation theorem. In general, a variety of different physical 
mechanisms can contribute to damping, including bulk and 
surface defects, thermoelastic damping, nonlinear coupling to 
other modes, phonon-electron coupling, clamping loss etc. 
Identifying experimentally the contributing mechanisms in a 
given system can be highly challenging, as the dependence 
on a variety of parameters has to be examined systematically. 
Nanomechanical systems suffer from low quality factors Q 
relative to their macroscopic counterparts [2]. This behavior 
suggests that damping in nanomechanical devices is dominated 
by surface properties, since the relative number of atoms 
on the surface or close to the surface increases as device 
dimensions decrease. This point of view is also supported by 
some experiments [9], [10]. However, very little is currently 

The authors are with the Electrical Engineering Department, Technion, 
Israel Institute of Technology. 



known about the underlying physical mechanisms contributing 
to damping in these devices. 

In the present paper we study damping in a nanomechanical 
oscillator operated in the nonlinear regime. Nonlinear effects 
are of great importance for nanomechanical devices. The 
relatively small applied forces needed for driving a nanome- 
chanical oscillator into the nonlinear regime is usually easily 
accessible. Thus, a variety of useful applications such as 
frequency synchronization, frequency mixing and conversion, 
parametric and intermodulation amplification [11], mechanical 
noise squeezing [12], and enhanced sensitivity mass detection 
[13] can be implemented by applying modest driving forces. 
Moreover, monitoring the displacement of a nanomechanical 
resonator oscillating in the linear regime may be difficult 
when a displacement detector with high sensitivity is not 
available. Thus, in many cases the nonlinear regime is the 
only useful regime of operation. However, to optimize the 
properties of NEMS devices operating in the nonlinear regime 
it is important to characterize the effect of damping in this 
regime. 

The effect of nonlinear damping for the case of strictly 
dissipative force, being proportional to the velocity to the p'th 
power, on the response and bifurcations of driven Duffing 
[14], [15], [16], [17] and other types of nonlinear oscilla- 
tors [18], [19], [16] have been studied extensively. For the 
present case we consider a Duffing oscillator having nonlinear 
damping force proportional to the velocity cubed. As will be 
shown below, this approach is equivalent to the case where 
damping nonlinearity proportional to the velocity multiplied 
by the displacement squared is considered (see Ref. [20]). 
We have recently studied a closely related problem of a 
nonlinear stripline superconducting electromagnetic oscillator 
[21], [22], where nonlinear damping was taken into account. 
With some adjustments, these results are implemented for the 
case of a nanomechanical nonlinear oscillator. To determine 
experimentally the rate of nonlinear damping, as well as the 
Kerr constant and other important parameters, we measure the 
response near the resonance in the nonlinear regime [23], [24]. 
Measuring these parameters under varying conditions provides 
important insights into the underlying physical mechanisms. 

II. Experimental setup 

For the experiments we employ nanomechanical oscillators 
in the form of doubly clamped beams made of PdAu (see 
Fig. The bulk nano-machining process used for sample 
fabrication is similar to the one described in Ref. [24]. 
The dimensions of the beams are length 100-200/im, width 
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0.25- 1/im and thickness 0.2/im, and the gap separating the 
beam and the electrode is 5^m. Measurements of mechanical 
properties are done in-situ a scanning electron microscope, 
where the imaging system of the microscope is employed 
for displacement detection [24]. Some of the samples were 
also measured using an optical displacement detection system 
described elsewhere [12]. Driving force is applied to the beam 
by applying a voltage to the nearby electrode. With a relatively 
modest driving force the system is driven into the region of 
nonlinear oscillations [24], [25]. 




Fig. 1. The device consists of a narrow cantilever beam (length 200/^m, 
width l-0.25^tm and thickness 0.2^m) and wide electrode. The excitation 
force is applied as voltage between the beam and the electrode. 



III. Equation of motion 

We excite the system close to its linear fundamental mode. 
Ignoring all higher modes allows us to describe the dynamics 
using a single degree of freedom x. 

The nonlinear equation of motion is 



mx + 2bii + k\X + b 3 x 3 + k 3 x 3 



cap 



dx 



where m is the effective mass of the fundamental mode, 
£cap = C(x)V 2 /2 is the capacitance energy, C{x) — Cq/{1 — 
x/d) is the displacement dependent capacitance, d is the gap 
between the electrode and the beam, b\ is the linear damping 
constant, 03 is the nonlinear damping constant, k\ is the linear 
spring constant and k 3 is the nonlinear (Kerr) spring constant. 

The applied voltage is composed of large DC and small AC 
components V(t) = Vbc +wcos(wt) where v is constant and 
v <C Vdc- Thus the equation of motion reads 

(1 - x/d) 2 (x + 271X + to 2 x + 7 3 i 3 + a 3 x 3 ) 

= F (1 + 2f cos(ut)) , (2) 



where u>q — y kifm, 71 = bi/m = ujo/2Q (Q being the 
mechanical quality factor), 73 = 63 /m, 0:3 — k 3 /m, F = 
CoVp C /2md and / = v/Vdc- 

IV. Multiple scales approximation 

We use the standard multiple scales method to solve Eq. 
[18], [26]. The harmonic excitation frequency is assumed to 
be close to the primary resonance 

U) = UJ + (7, 



where a <C luq is a small detuning parameter. We also 
assume the linear damping coefficient 71, both coefficients of 
nonlinear terms, 73, 0:3, and 1 jd to be small. Keeping terms up 
to the first order in the small parameters leads to the following 
form of the solution for x{t) 



F 



x{t) = — + {A (t) e^ ot + c.c.) , 



(3) 



where the first term in the right-hand side is a constant 
displacement due to electrostatic attractive force, and A (t) 
is slowly varying envelope (on the time scale of 1/ojq). The 
differential equation for A (i) in this approximation is given 

by 



2U! 



dA ,. . . . 
1— + («7i + Aw ) A 



dt 



+ 3 (a 3 + H3Wo) A 2 A* = Ffe 1 



(4) 



where A^o = (SasF/uiQ — 2/d) F/2ojq is a small correction 
to the linear resonance frequency ujq. 

The solution for A (t) can be represented as 



A(t) ^ae l ^ +AuJot \ 



(5) 



where a and <p are rea l- Substituting Eq. (|5} into Eq. (0} and 
separating real and imaginary parts one finds 

-2cj a— + 3a 3 a 3 = FfcosU - Atut) , (6a) 
dt 



— n ( ^ + 7i a 



3731^ = Ff sin (0 - Aujt) , (6b) 



where Ao; = a — Acjq is the excitation frequency detuning 
from the shifted resonance frequency loq + Aujq. In the steady 
state a and <fi — Aujt are constant and the following equation 
for the steady state response amplitude a can be derived from 
Eq. © 

9 (aj + 7 3 2 c^) a 6 + 12uj (7173^0 - A ^a 3 ) a 4 

+ Auo 2 (Auj 2 + 7l 2 ) a 2 - F 2 f 2 = 0, (7) 

Equation of the same form was obtained in Ref. [21], 
where a superconducting oscillator having Kerr nonlinearity in 
addition to nonlinear damping was considered. All subsequent 
analysis is thus based on Ref. [21]. 

When 73 is sufficiently small the solutions of Eq. |7} behave 
very much like the ordinary Duffing equation solutions to 
which Eq. Q reduces to when 63 = (see Fig. |2j- 

Interestingly enough, equations similar to Eq. @ and Eq. 
arise when the damping nonlinearity is considered to 
be proportional to velocity multiplied by the displacement 
squared (instead of velocity cubed) 



mx + 2b\x + k\x + 63X x + k 3 x 



d£ c 



dx 



(71 +73^ 2 / 2w o) and 73 by 
Eq. (0J and Eq. Q gives the 



Substituting ji by 
73 (1 + IF/lu 2 ) /3lu 2 i 
correct relations for this case. Therefore, the behavior for 
these two cases is similar near the resonance frequency. 
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Fig. 2. Steady state solutions under different excitation amplitudes /. In case 
f < fc only one real solution exists, no bistability is possible. In case / = f c 
the system is on the edge of bistability, and one point exists where a 2 vs. ui 
has infinite slope. In case / > f c the system is in bistable regime having 
three real solutions over some range of frequencies. Two of these solutions 
are stable. 



V. Special points 

Referring to Fig. 0we define some points in the a 2 vs. co 
curves which we use in experimental data analysis. 

The first point is the maximum response, shifted by Aw m 
from u>o + Auiq an d having the amplitude a m . Differentiating 
Eq. Q with respect to Aui and demanding d (a 2 ) /dAaj = 
yields 



2uiqAu) ti 
3a 3 



(8) 



Another point of special interest is the point where 
the jump in amplitude occurs and therefore the condition 
dAw/d (a 2 ) = must be satisfied. Applying this condition 
to Eq. Q yields 



27 (a 2 + t£w§) a 4 + 2Au a (7i73^o - AuJa s) a 2 



(Auj 2 



7?) 



= 0. 



■j ^ -r-rij — (9) 

Eq. has a single real solution at the point of critical 
frequency Auj c and critical amplitude a c , where the system 
is on the edge of bistability. This point is defined by two 
conditions 

dAuj 



d(o 2 ) 
d 2 Aw 



= 0, 



Q. 



d(a 2 ) 2 

In general, 73 is positive but a 3 can be either positive (hard 
spring) or negative (soft spring). In our experiment a 3 > 0. 
By applying these conditions one finds 



Aw, 



71 P + 3 
7i^o 1 



3V3 "3 l-p' 



(Ha) 
(lib) 



where p = \/373Wo/q;3. The driving force at this critical 
point is denoted in Fig. [2] as fc- Note that bistable region 
is accessible only when p < 1. 

VI. Experimental data and results 

A typical measured response of the fundamental mode of 
a 200/im (125/im) long beam occurring at /o = 123.2 kHz 
(/o = 524.6 kHz) measured with Vdc = 20 V and varying 
excitation amplitude is seen in Fig. 0a) (Fig. 0b)). We derive 
the value of 71 = uiq/2Q from the linear response at low 
excitation amplitude and find Q = 7200 (Q = 2100). 
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Fig. 3. Measured response vs. frequency shown for both upward and 
downward frequency sweeps with Vdc = 20V and with varying peak-to- 
peak excitation amplitude V pp . (a) 200/xm long beam with fundamental mode 
occuring at fo = 123.2 kHz. (b) 125^tm long beam with fundamental mode 
occuring at /q = 524.6 kHz. 



Our displacement detector is highly nonlinear, introducing 
thus a significant distortion in the measured response. In 
order to minimize the resultant inaccuracies, we employ the 
following method to extract the nonlinear parameters. 

In general, the sum of the three solutions for a 2 at any 
given frequency can be found from Eq. @. This is employed 
for the jump point at ujq + Aujj seen in Fig. Using Eq. (0 
to calibrate the measured response at this jump point one has 



{2h 1 + h 2 ) 



2uiQAuj ri 
3a.3 



Aujq (71 73^0 - Acjj-a 3 ) 



3 (a 2 



- 7 3 2 ^o 6 ) 
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Fig. 4. Experimental results forp = t/373^o / a 3 vs - peak-to-peak excitation 
amplitude V pp . (a) 200^tm long beam with fundamental mode occuring at 
/o = 123.2 kHz and Q = 7200. (b) 125/im long beam with fundamental 
mode occuring at /o = 524.6 kHz and Q = 2100. 



or 



(2hi + h 2 ) Aw r , 



1 



71 



V3 



Au 



= 0, (13) 



where hi and hi are defined in Fig. [2] Due to the frequency 
proximity between the maximum point and the jump point at 
OJ = luq + Auij the inaccuracy of such a calibration is small. 
Moreover, as long as excitation amplitude is high enough, h 2 is 
much smaller than hi and even considerable inaccuracy in h 2 
estimation will not have any significant impact. This equation 
can be used to estimate p for different excitation amplitudes. 
The results of applying Eq. Jl 3i to experimental data from the 
two different beams can be seen in Fig. |4] 

Using this procedure we find p « 0.35 for the 200/im 
long beam and p as 0.14 for the 125/im long beam. These 
results are identical to the values of p estimated using Eq. 
i ll lab . Referring to Eq. Jl li and Eq. l|9) we see that in our 
system the damping nonlinearity is not negligible and has a 
measurable impact on both the amplitude and frequency offset 
of the critical point, as well as on jump points in the bistable 
region. 

To determine the value of 0:3 we measure the static de- 
flection of the beam's center as a function of an applied 
DC voltage Vdc- From a fit to theory we find a.3 = uiq ■ 
0.092 /1m" 2 for the 200/im long beam. 

VII. Conclusions 

In this work we have demonstrated conclusively that non- 
linear damping in nanomechanical doubly-clamped beam os- 
cillator may play an important role. The method presented in 
this paper may allow a systematic study of nonlinear damping 
in nano-mechanical oscillators, which may help revealing the 
underlying physical mechanisms. 
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